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ABSTRACT 



Aims. We study the variability of the Fe 6.4 KeV emission line from the Class I young stellar object Elias 29 in the p Oph cloud. 
Methods. We analysed the data from Elias 29 collected by XMM-Newton during a nine-day, nearly continuous observation of the 
p Oph star-forming region (the Deep Rho-Oph X-ray Observation, named Droxo). The data were subdivided into six homogeneous 
time intervals, and the six resulting spectra were individually analysed 

Results. We detect significant variability in the equivalent width of the Fe 6.4 keV emission line from Elias 29. The 6.4 keV line 
is absent during the first time interval of observation and appears at its maximum strength during the second time interval (90 ks 
after Elias 29 undergoes a strong flare). The X-ray thermal emission is unchanged between the two observation segments, while 
line variability is present at a 99.9% confidence level. Given the significant line variability in the absence of variations in the X-ray 
ionising continuum and the weakness of the photoionising continuum from the star's thermal X-ray emission, we suggest that the 
fluorescence may be induced by collisional ionisation from an (unseen) population of non-thermal electrons. We speculate on the 
possibility that the electrons are accelerated in a reconnection event of a magnetically confined accretion loop, connecting the young 
star to its circumstellar disk. 
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1. Introduction 

The X-ray emission from young stellar objects (YSOs) at CCD- 
resolution is usually modelled as thermal emission from a hot 
plasma in coronal equilibrium, with higher characteristic tem- 
peratures than observed in older and less active stars. An inter- 
esting deviation from a pure thermal X-ray spectrum is the pres- 
ence of fluorescent emission from neutral (or weakly ionised) 
Fe as show n by the presen c e of th e 6.4 keV line. This was first 
detected by Imanis hi et al.l d2001 .) in the X-ray emission of the 
YSO YLW16A in p-Oph, during a large flare: in addition to the 
Fe xxv complex at 6.7 keV, a 6.4 keV emission line was clearly 
visibe. Such fluorescence line is produced when energetic X- 
rays photoionise cold material close to the X-ray source, and it 
is therefore a useful diagnostic tool of the geometry of the X-ray 
emitting source and its surroundings. 

Since 2001, detections of the Fe K fluorescent emission line 
at 6.4-keV in the spectra of YSO s have been reported by a 
number of authors. iTsuiimoto et al.l d2005l) has identified seven 
sources with an excess emission at 6.4 keV among 127 ob- 
servations of YSOs within the COUP observation of Orion; 
Fava ta et"al1 d2005bl) report 6.4 keV fluorescent emission in Elias 
29 in p-Oph both during quiescent and flaring emission, unlike 
all other reported detection of Fe fluores cent emission in YSOs 
that were made during intense flaring; iGiardino et alj d2007l) 
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have detected Fe 6.4 keV emission from a low-mass young 
star in Serpens, during an intense, long-duration flare. Recently, 
ICzesla & Schmittl d2007l) have reported intense Fe fluorescent 
emission in the spectrum of V 1486 Ori during a strong flare, 
when the plasma reached a temperature in excess of 10 keV. 

The 6.4 keV fluorescent line has been detected in different 
classes of X-ray emitters: X-ray binaries, active galactic nuclei 
(AGNs), massive stars, supernova remnants, and the Sun itself 
during flares. In the case of the Sun, the fluorescing material is 
the solar photosphere, in the YSOs, however, indications are that 
the material in the circumstellar disk and its related accretion 
structures could be responsible for the fluorescence. The typi- 
cal equivalent width of the 6.4 keV emission line in the stud- 
ies mentioned above is of the order of 150 eV and is too large 
to be explained with fluorescent emission in t he stellar photo- 
sphere or in diffuse circu mstellar material (e.g. Tsuiimot o et all 
12001 iFavata et all l2005bTl . This scenario implies that the disk 
is "bathed" in high-energy X-rays emitted by the star, with sig- 
nificant astrophysical implications; for instance, X-rays, in addi- 
tion to cosmic rays, would play an important role in photoion- 
ising the circumstellar material around young star and thus in 
coupl ing the gas to the am b ient magnetic field (a s suggested by 
e.g. lGlassgold et alll2000l) . ICeccarelli et alj d2002l) suggest that 
the "hot" component they observe in the disk, in the infrared, is 
heated by the stellar high-energy emission. 
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Fig. 1. The light curve of Elias 29 over the nine days of the Droxo observation from PN (top), MOS1 (bottom-left), and MOS2 
(bottom-right). The line with error bars gives the background-subtracted light curve of the source, while the thinner line without 
error bars gives the total counts (source plus background). The 6 time intervals that we selected for the spectral analysis, on the 
basis of the PN data, are indicated. 



Observations of the time variability of the Fe 6.4 keV emis- 
sion in YSOs would provide useful constraints on the geometry 
and sizes of the star-accretion disk system and of other circum- 
stellar structures such as funnel flows, jets, or wind columns. We 
present here the results of a time-resolved spectral study of the 
X-ray emission of Elias 29, during ~ 9 days of nearly contin- 
uous observation by XMM-Newton in the context of the ultra- 
deep observation of p-Oph, named Droxo (from Deep Rho-Oph 
X-ray observation - iPillitteri et all l2007t) . We investigated the 
presence of variations in its strong Fe 6.4 keV emission over the 
9-day time scale covered by the observations. 

This paper is structured as follows. After a summary, below, 
of the properties of Elias 29, the observations and data analy- 



sis are briefly presented in Sect. [2] Results are summarised in 
Sect.|3j the simulations carried out to assess the reliability of 
the line detections and the significance of its variability are de- 
scribed in Sect. [4] The results are discussed in Sect. [5] 



1.1. Elias 29 (GY214) 

Elias 29 (16:27:09.4, -24: 37:18.9) with a bolome t ric luminos- 
ity L = 26 - 27.5 Lq dBontemps et al.L 1200 U iNatta et all 
2006 ) is the most lum inous Class I YSO in the p-Oph cloud. 
Muz erolle et al.l dl998l) used the luminosity in the Bry line to 
determine the object's accretion luminosity at L acc = 15 - 18L0, 
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which makes it the source wit h the highes t accre tion luminosity 
in their sample. More recently. iNatta et ail (12006) used the lumi- 
nosity of the hydrogen recombination lines to derive an accretion 
luminosity of 28.8 Lq. 

Us ing millimeter interferometric observations, Boogert et al.l 
(2002) resolved the emission from the disk and the envelope sur- 
rounding Elias 29, showing that the disk is in a relatively face- 
on orientation (i < 60°), which explains many of the remarkable 
observational features of this source, such as its fiat spectral en- 
ergy distribution, its brightness in the near-infrared, the extended 
components found in speckle interferometry observations, and 
its high-velocity molecular outflow. Their best-fitting disk model 
has an inner radius of 0.01 AU, outer radius of 500 AU, and a 
mass M = 0.012M©. 

Elias 29 was previously observed in X-rays with ASCA, 
Chandra and XM M-Newton. In the Chandra observation 
dlmanishi et al.L 1200 ll) . the source quiescent phase is charac- 
terised by a temperature of 4.3 keV and luminosity of 2.0 x 
10 30 erg s _1 , fully consistent with the val ues derived from th e 
subsequent XMM-Newton observations by Oza wa etaf] (|2005): 
kT = (3.6 - 5.1) keV, N(H) = (4.4 - 5.3) x 10 22 cm" 2 , Z = 
(0.8 - 1.3) Zq, and Lx = 2.8 x 10 30 erg s _1 . The s ource was seen 
flarin g during one of the ASCA observations (Kamata et al., 
1 1997b and during the Chandra observation. The two flares had 
similar intensity and duration with an e-fo lding time of ~ 10 ks 
dTsuboi et al.L 120001: llmanishi et al.U2001l) . 

2. Observations and data analysis 

The Droxo program is a nominal 500 ks observation of the 
p-Ophiuchi star-forming region performed by the EPIC cam- 
era on board the XMM-Newton satellite. The observation was 
performed over 9.4 days, starting 8 March 2005 (orbits 0961- 
0965). Detai ls of the observations and data reduction procedure 
are given in iPillitteri et al.l d2007). The preliminary data reduc- 
tion was done with SAS software version 6.5 in order to obtain 
lists of photon events calibrated both in energy and astrometry 
for the three instruments, MOS1, MOS2, and PN, for each or- 
bit. The data were filtered in the energy band 0.3-10. keV, and 
only events with pattern < 4 and flag = were retained for the 
spectral analysis. The spectral analysis was performed using the 
xspec package VI 1.2, after rebinning the spectra to a minimum 
of 20 source counts per (variable width) spectral bin. 

3. Results 

Figure Q] shows the light curve of Elias 29 from the three instru- 
ments, PN, MOS1, and MOS2, and the six time intervals that 
we selected for our spectral analysis. On the basis of the PN 
data, we selected 5 time intervals with low background (here- 
after "segl" to "seg5"), plus one time interval covering the 
strong flare at about 94 ks from the beginning of the observa- 
tion. During "segl", the MOS1 camera was very likely hit by a 
micro-meteorite (which compromised one of the chips) and the 
instrument was switched off for the rest of this segment of the 
observation, during the flare and part of "seg2". MOS2 data are 
available for all the segments of observation, but are insufficient 
during the flare, thus, we did not include the "flare" time interval 
in the joint spectral analysis of PN and MOS data. In addition, 
we only used PN and MOS2 data for "segl" and "seg2". 

The spectra were initially modelled by an absorbed one- 
temperature plasma model. The results of the six spectral fits 
for the PN data alone are reported in Table Q] and the results of 



the (five) joint fits to the PN, MOS1, and MOS2 data in Table 
[2] For each time-interval, the values of the spectral parameters 
derived from the simultaneous fitting are very similar to the val- 
ues derived from the PN data alone. The error bars in the best-fit 
parameters improve marginally, but^ 2 values worsen, likely be- 
cause of calibration uncertainties. 

The spectra and the fits for "segl" and "seg2" are shown 
in Fig. [2] and the spectrum for the flare in Fig.[3] The aver- 
age values of the source's spectral parameters while quiescent 
(N(H) = 6.8 x 10 22 cm- 2 , kT = 3.7 keV, and Z = 0.8 Z©) are 
very similar to the values derived from previous observations, so 
is its quiescent luminosity Lx ~ 10 30 erg s _1 , with no evidence 
of long-term variability. 

The X-ray flare is similar to other events previously observed 
from this source. During the flare, the source counts first in- 
creased impulsively by a factor of ~ 8 and then decreased ex- 
ponentially with a decay time of ~ 6 ks. The source spectrum 
during the flare is shown in Fig. [3] together with the spectral fit. 
As presented in TableQ] the fitted value of the plasma temper- 
ature does not appear to change significantly during the flare, 
possibly due to the stringent processing criteria applied to the 
data that resulted in the events of the flare peak being discarded. 

From Tables[T]and|2] it is apparent that the best-fit values of 
N(H), kT, and Z for the different time intervals do not show sig- 
nificant variations, since they are all consistent with each other 
within 2<x. This is also apparent by comparing the PN and MOS2 
spectra from "segl" and "seg2" in Fig. [2] where the overall spec- 
tral shape and amplitude are very similar during the two time 
intervals (in both instruments). On closer inspection, however, 
a significant difference between the two time intervals becomes 
apparent: during time interval "seg2", a visible excess of emis- 
sion around 6.4 keV, the energy of the Fe fluorescent line, is 
present both in the PN and MOS2 data. 

To quantify this excess and monitor its variation, we repeated 
the spectral fits of the spectra with an absorbed IT plasma model 
and an additional Gaussian line component at 6.4 keV. The posi- 
tion and width (10 eVQ) of this component were constrained dur- 
ing the fit, while its normalisation was left free to vary. The other 
parameters of the absorbed IT plasma model (absorbing column 
density, temperature, and normalisation of the thermal spectrum) 
were also free parameters in the fit. The spectra and spectral fits 
to the six PN spectra in the energy range AE =4-8 keV are 
shown in Fig. [4] The strong excess emission at 6.4 keV in the 
spectra of "seg2" is well accounted for by the fitted line, while 
it is clear that no such additional line at 6.4 keV is needed to 
fit the data from "segl". The four other spectra are somewhat in 
between these two extremes in regard to an excess of emission 
at 6.4 keV. 

The results from the spectral fits to the PN data with the ad- 
ditional line at 6.4 keV are summarised in Table [3] As can be 
seen by comparing Tables Q] and [3] the fitted values for absorb- 
ing column density, temperature, metallicity, and normalisation 
of the thermal spectrum are not affected by the addition of the 6.4 
keV-line component. The null-hypothesis probability of all the 
fits improves, except for "segl", for which the null-hypothesis 
probability decreases marginally. Indeed, an excess of emission 
at 6.4 keV does not appear to be present at all in this spectrum, 
so that there is no reason to add an extra line component at 6.4 
keV. We have also added this component to the spectral fit of 
"segl" to make the fits homogeneous and directly comparable. 



1 The iron Ko- fluorescence line consists of two components with an 
energy separation of ~ 10 eV: Ko-i and K«2 at 6.404 and 6.391 keV, 
respectively, for Fe i 
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Fig. 3. The spectrum of Elias 29 (PN data) during the strong 
flare between "segl" and "seg2", together with the best-fit one- 
temperature model. 



The equivalent width of the fitted line appears to vary signifi- 
cantly from the 13 eV of "segl" (where the excess emission at 
6.4 keV does not seem to be present) to the 249 eV of "seg2". 

Table|4] summarises the results of the simultaneous spectral 
fitting of PN and MOS data with an absorbed IT plasma model 
and an additional Gaussian line component at 6.4 keV: the results 
obtained from this joint spectral analysis are again fully consis- 
tent with the results from the PN data alone (compare Tableland 
@). The strong variation in the intensity of the excess emission 
at 6.4 keV between "segl" and "seg2" is confirmed; the best-fit 
line equivalent width varies from 10 eV in "segl" (where the line 
also appears to be absent in the MOS2 data) to 194 eV in "seg2" 
(where the excess of emission is also strong in the MOS2 data 
- see Fig.|2]i. For "seg3" and "seg4", the equivalent width of the 
emission at 6.4 keV derived from the joint fits differs by ~50% 
from the value derived from the PN data alone, but best-fit val- 
ues of the line's flux are consistent within the two data-sets; for 
"seg5", the equivalent width of the emission at 6.4 keV derived 
from the two spectral analysis are very similar. 

As described in Sect.|4] we have performed a set of Monte 
Carlo simulations of the PN spectra to assess whether the line at 
6.4 keV seen in the spectra from "seg2" to "seg5" is real rather 
than the result of statistical fluctuations in the data. The results of 
these simulations are summarised in the last column of Table|3] 
which gives the probability that a line at 6.4 keV, with equivalent 
width equal to, or greater than, the best-fit value, could be the re- 
sult of random fluctuations for each segment. This probability is 
high for "segl" and the "flare" segment, for which the data are 
consistent with the absence of a line in the spectrum. Conversely, 
the probability that the excess of emission in "seg2" and "seg5" 
is due to noise fluctuations is low (0.1% and 0.2%, respectively). 
For time intervals "seg3" and "seg4" the situation is less clear, 
since with a probability of 4.2% and 3.5%, respectively, that the 
excess at 6.4 keV could be due to random fluctuations, the evi- 
dence for intrinsic emission is less compelling. Nevertheless an 
analysis of the spectrum of the source integrated over the time 
intervals "seg3", "seg4", and "seg5" shows that the evidence for 
the presence of the line is very strong, with a probability that the 
excess at 6.4 keV is due to random fluctuation of less than one 
in a thousand (as summarised in the last line of Table[3]l. 

The simulations described in Sect.|4]also show that the vari- 
ation in the line's equivalent width seen between "segl" (where 




channel energy (keV) 

Fig. 5. The spectrum (PN data) of Elias 29, between 4 and 8 
keV, integrated over the four time intervals of quiescent emis- 
sion where an excess at 6.4 keV appears to be present. With the 
improved statistic of this long time span, the presence of 6.4 keV 
Fe fluorescent emission is very clear. 



the line appear to be absent or extremely weak) and "seg2" are 
very likely intrinsic (with a probability of 99.9%). 

To further constrain the properties of the Fe 6.4 keV emis- 
sion from Elias 29, we derived the spectrum of the source inte- 
grated over the quiescent time intervals where the excess at 6.4 
keV appears to be present (although with varying significance), 
i.e. intervals "seg2", "seg3", "seg4", and "seg5" (PN data only). 
The summed spectrum between 4 and 8 keV is shown in Fig. [5] 
where the emission line at 6.4 keV is very clear thanks to the 
the improved statistic. A fit with an absorbed IT model to this 
spectrum yields N(H) = 6.7 x 10 22 cnr 2 , kT = 3.7 keV and 
Z = 0.7 Zq with a null hypothesis probability P — 3.5 x 10~ 5 . 
The spectral parameters are very similar to the average values 
derived above from the fits to the individual five quiescent time 
intervals, although the fit probability is rather low. With the ad- 
dition of a line at 6.4 keV, the fit probability increases substan- 
tially to P — 2.1 x 10~ 3 (and the values of the parameters of the 
IT model hardly change). In this case we left the position of the 
Gaussian line free to vary during the fit, together with its normal- 
isation. The fitted line equivalent width is W(,a keV = 143 eV and 
its position is 6.44 + 0.03, fully consistent with being Fe K fluo- 
rescent emission. Indeed, although the reliability is limited (the 
quoted error is at lcr), the line peak energy of 6.44 keV would 
suggest that the fluorescing Fe is Neon-like: Fe xvn or so. The Fe 
K fluorescent emission line energy is, in fact, a slowly increasing 
function of ion isation state, rising from 6.40 keV in F e i to 6.45 
keV in Fe xvn (iHousel [19691 iGeoree & Fab"ianL[T99l . 

4. Simulations 

To quantify the reliability of the detection of the 6.4 keV emis- 
sion in the different spectra of Elias 29, we performed a set of 
Monte Carlo simulations of the PN spectra. For each time in- 
terval discussed in Sect. [3] 1000 random realisations of the ab- 
sorbed IT- model fitted to the spectrum without the additional 
Gaussian component at 6.4 keV (as per Table[TJ were gener- 
ated using xspec These random realisations have the same noise 
characteristics of the data and the same energy binning of the 
real data (for which channels are combined to a minimum of 20 
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Fig. 2. Elias 29 spectra from PN (fop) and MOS2 (bottom) with spectral fits, during "segl" (before the flare, left) and "seg2" (after 
the flare, right). The spectra (from the two time intervals) are very similar in overall shape and amplitude, and the fits (with an 
absorbed IT plasma model) result in very similar parameters. During "seg2", however, a significant excess of emission at 6.4 keV 
is present, which is not visible in the data from "segl". All the spectra are rebinned to a minimum of 20 source counts per (variable 
width) spectral bin. 



counts per energy bin). The 1000 simulated spectra were then fit- 
ted with the same procedure used for the real data. First a fit with 
an absorbed IT plasma model was performed then a Gaussian 
line component at 6.4 keV with a narrow cr of 10 eV was added 
to the model and the simulated data were re-fitted. The line po- 
sition and width were kept fixed during this fit, while its normal- 
isation was left unconstrained. The equivalent width of the fitted 
line was then computed for all 1000 fits for each time interval. 

Since the input model to the simulated spectra did not con- 
tain any emission line at 6.4 keV, the number (over 1000) of fit- 
ted spectra with a fitted line with a given W(,a keV gives the prob- 
ability that a line with that equivalent width could be measured 
in the real spectrum because of the noise's random fluctuations. 
For example, in the case of "seg 1" the number of simulated spec- 
tra with W6.4kev ^ 13 eV (W6.4kev = 13 eV being the best-fit 
equivalent width in the real data) is 453, thus, the probability of 
deriving a spurious line with W(,a keV ^ 13 eV from the spec- 
trum of "segl" is 45%. This can be seen by looking at the top 
panel of Fig. [6] which gives the distribution of W(,a keV for the 
fitted lines for the simulation of "segl". Note that, in the case 
of "segl", there is no visible excess of emission at 6.4 keV, so 
a high incidence of simulated spectra with W(,a keV ^ 13 eV is 
consistent with the absence of the line in the intrinsic spectrum. 

In the case of "seg2", on the other hand, the number of sim- 
ulated spectra with W(,a keV ^ 249 eV is 1 (over 1000 simula- 



tions), indicating that the excess at 6.4 keV seen in the spectrum 
of "seg2" is very likely intrinsic to the source. The results of 
all these simulations for all the different time intervals, in terms 
of the probability of by chance observing a line with a value 
of W(,a keV greater or equal to the one measured in the data, are 
summarised in Table|3] For illustration, the distribution of the fit- 
ted spectral parameters N(H), kT, and Z for the 1000 simulations 
for "segl" are given in Fig. [7] 



Finally, to assess the probability that a line with the same 
equivalent width as derived in "seg2" (WI5.4 keV = 249 eV) could 
be present in the source's intrinsic spectrum of "segl", but could 
go undetected due to the noise fluctuation, we performed one 
additional simulation for "segl". In this case, the simulation in- 
put model was the best-fit absorbed IT model for "segl" (as per 
Table[TJ with the addition of a Gaussian line at 6.4 keV with 
W(,a keV = 249 eV. In this case, the number of simulated spectra 
for which we derived a fitted line at 6.4 keV with W(,a keV ^ 13 
eV is 1 (over 1000), as can be seen from Fig. [6] In other words, 
the probability that during "segl" an intrinsic emission line with 
the same W(,a keV as derived in "seg2" could result in an observed 
excess with W(,a keV ^ 13 is 0.1%. This indicates that the strong 
variations in Fe 6.4 keV emission detected in the data between 
"segl" and "seg2" are very likely intrinsic. 




5 ^5 

channel energy (keV) channel energy (keV) 

seg2 seg3 
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Fig. 4. Spectra (PN) and spectral fits of Elias 29 during 6 different time intervals, between 4 and 8 ke V. From left to right, top to 
bottom, spectra are from time intervals "segl", "flare", "seg2", "seg3", "seg4", and "seg5" (these intervals are indicated in Fig. 
1). Fits were performed using an absorbed IT thermal model with the addition of a line at 6.4 keV (indicated in the plots by a 
dash-dotted line). Note the strong excess of emission at 6.4 keV in the spectra for "seg2", which does not appear to be present in the 
spectrum from the first time interval. 



5. Discussion 

As mentioned in the introduction, previous works have consid- 
ered that the most likely location of the cold material traced by 
fluorescent emission is, in YSOs, the circumstellar disk, as no 
fluorescence, whether from the photosphere or from a circum- 



stellar shell, can result in the large equivalent widths observed 
(approximately 150 eV), given a plausible thermal photoionis- 
ing continuum. 

The radiative transfer computation necessary for determin- 
ing the fluorescent equivalent width in the case of an optically 
thick target (i.e. a stellar photosphere) shows that, in the case 



G. Giardino et al.: Results from Droxo. I, The variability of fluorescent Fe 6.4 keV emission in the young star Elias 29 



7 




60 




Fig. 6. Distributions of the best-fit equivalent width of the 6.4 
keV line in two sets of simulations of the spectra of "segl". 
In the top panel the input model to the simulation did not con- 
tain a Gaussian line component at 6.4 keV, so the distribution 
gives the probability of deriving a spurious emission line with a 
given equivalent width. The best-fit equivalent width derived for 
"segl" of 13 eV is consistent with the line being spurious. In the 
bottom panel the input model contained a Gaussian line compo- 
nent at 6.4 keV with W(,a k e v = 249 eV, as derived in the data 
for "seg2". In this case, the probability of detecting a Gaussian 
component with W(,a k e v ^ 13 eV is very low (0.1%). Each sim- 
ulation consisted of 1000 random realisations of the input model. 



of a semi-isotropically illuminated slab (such as a photosphere), 
the equivalent width W6.4kev varies approximately from 150 
to 90 eV (assuming photospheric solar iron abundance), for a 
power- law inciden t spect rum with spectral index T = 1 .2 to 2. 1 
dGeorge & Fabiani 119911) . Given that the thermal spectrum of 
stellar sources is g enerally steeper than a power-law spectrum 
with T ~ 2.0 (e.g. iFavata et all l2005bt iGiardino et alt 120071) . 
photospheric emission can hardly explain the observed equiva- 
len t widths. 

iGeorge & Fabiani dl99ll) show that higher equivalent widths 
for the fluorescent Fe K line can be obtained if the cold mate- 
rial is distributed in a centrally illuminated disk. In this case, for 
favourable disk inclination angles (near face-on, i = 0° - 30°), 
W6.4kev varies approximately between 120 and 200 eV for an 
incident spectrum spectral index F = 2.3 to 1.3. This interpreta- 
tion of th e origin of the fluore scent Fe K emission was also pro- 
posed bv lFavata et al.l ((2005b) for Elias 29, whose accretion disk 
was shown to be in a relatively face-on orientation (i < 60°) by 
Boog ert et alj d2002l) . This interpretation cannot, however, easily 
explain the current set of data. 

The remarkable feature of the present data is that the strong 
Fe 6.4 keV emission, apparent in the second segment of the ob- 




kT [keV] 




Fig. 7. Distributions of the best-fit values of the absorbing col- 
umn density, plasma temperature, and metallicity in the the sim- 
ulations of the spectra of "segl". The vertical dashed lines give 
the parameter values used in the simulation's input model. The 
simulation consisted in 1000 random realisations of the spectral 
model fitted to the data of "segl". 



servation, is absent (or significantly weaker) during the first seg- 
ment of the observation, while the observed thermal continuum 
X-ray spectrum is basically unchanged. Although a stellar flare 
takes places between the two intervals, the radiative lifetimes of 
fluorescent Fe transitions are negligible compared to the time be- 
twe en the flare and the se cond observation segment (91 ks) - see 
e.g. Langhans et al.1 (1 19951) for the radiative lifetimes of Fe i and 
Schn abel et alj ( 2004 ) for Fe n - and therefore the flare cannot 
explain the 6.4 keV emission observed in the second segment of 
the observation. To explain the observation, a sustained mecha- 
nism ionising the neutral Fe atoms must operate during the rest 
of the observation, on time scale of days, independent of (and not 
correlated with) the thermal X-ray emission. This mechanism 
must be absent during the first segment. This is unlikely to be 
the stellar X-ray continuum emission, as this emission appears 
to be unchanged between the first and the second observation 
segments. 

A scenario in which a localised region of photoionised Fe 6.4 
keV emission is somehow eclipsed during "segl" is also difficult 
to reconcile with the data, given the high equivalent width of the 
line (W6.4kev ^ 150 eV), when this is present. As mentioned 
above, to explain such a large equivalent width by photoionisa- 
tion, the favorable geometry of a nearly face-on, centrally illumi- 
nated disk is needed in order to maximise the number of atoms 
whose line emission is free to escape. Requiring that the Fe 6.4 
keV emitting region is localised (and in order to be eclipsed by 
the star in a nearly face-on disk this volume needs to be signif- 
icantly localised), automatically reduces the equivalent widths 
that can be obtained. The same line of reasoning applies to a re- 
gion of localised emission on the stellar photosphere, even more 
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so, given that for the same incident spectrum one obtains smaller 
line equivalent widths. 

The large Fe 6.4 keV equivalent width derived in the sec- 
ond segment, 190 - 250 eV, is also difficult to reconcile with 
emission from photoioni sation of the mat e rial in the circumstel- 
lar disk. The model of George & Fabian] dl991l) for a centrally 
illuminated face-on disk and an incident X-ray spectrum with a 
power law index of F = 2.3, which is appropriate for Elias 2S0, 
predicts a maximum equivalent width of only 120 eV, that is, 70 
- 130 eV less than the observed best-fit value in segment 2. In 
addition, the maximum equivalent width of the 6.4 keV line cor- 
responds to the lowest temperature in the thermal spectrum, in 
contradiction with the expectations from a photoionisation pro- 
cess. 

An alternative explanation for the production of the Fe 6.4 
keV emission line in solar flares is collisional io nisation of K- 
shell electrons by a beam of non-thermal electrons dEmslie et al.L 
1986 ). This mechanism may also operate in some stellar sys- 
tems. lOstetTetal] d2007) have recently observed a superflare on 
the active binary system II Peg with the Swift telescope. Analysis 
of the X-ray spectrum from 0.8 to 200 keV from the XRT and 
BAT instruments reveals evidence of a thermal component with 
kT in excess of 7 keV, Fe 6.4 keV emission, and a tail of emis- 
sion out to 200 keV, which can be fit with a power-law model 
with spectral index r ~ 2. They attribute this power-law hard X- 
ray emission to non-thermal thick-target bremsstrahlung emis- 
sion from a population of accelerated electrons and find that 
collisional ionisation (of photospheric material) from the same 
non-thermal electrons is a more likely explanation for the 6.4 
keV feature than the photoionisation mechanism. 

In our case, a superflare is not taking place; indeed, no stellar 
flare is taking place during the second segment of observation. In 
a magnetospheric accretion scenario, however, ionised material 
is magnetically channeled in accretion stre ams that connect the 
star to the circumstellar accretion disk (e.g. lCalvet & Hartmannl 
119921) . The observation of very intense long-duration flares in 
YSOs, which implie s flaring loop lengths of tens of stellar radii 
dFavata et all 12005 ah . provides evidence of flaring associated to 
these accretion streams. 

A magnetic reconnection event leading to the acceleration 
of electrons within one of these structures, connecting the star 
to the disk, may or may not result in a (soft) X-ray flare event 
depending on the density of the material confined in the structure 
itself. If the density of material within the magnetic flux tube 
is not very high (with column densities from the acceleration 
site to the photosphere or to the disk of N £ 10 20 cm -2 ), the 
accelerated particles will impact the stellar photosphere and heat 
and "evaporate" the chromospheric plasma, providing the hot 
material that emits the X-ray flare observed between 0.1 and 10 
ke V. Indeed, th e flares f or which long flaring loops were derived 
in iFavata et al.l (12005 al) appear to take place in CTTs with no 
significant ongoing accretion (hence their "accretion" tubes are 
not mass-loaded). 

In YSOs, which are still strongly accreting, the density of 
material confined in the accretion tubes is expected to be high. 
For a fiducial case o f mass accretion rate of 1O" 7 M yr" 1 , 
iMuzerolle et all d200ll) estimate densities of n H = 10 12 - 10 13 
cm" 3 at ~ 2R„, that is, column densities of the order of N ~ 
10 23 - 10 24 cm" 2 assuming a stellar radius of 10 11 cm" 2 . In Elias 
29 with an accretion luminosity of 28.8 Lq, the a ccretion rate 
is estimated to be one order of magnitude higher dNatta et all 



2 as verified by performing a fit to the "seg2" spectrum of Elias 29 
with an absorbed power-law model 



2006), so the column density within the accretion tubes will 
be correspondingly higher. With such a high density, the acce- 
lareted electrons produced by a reconnection event in the ac- 
cretion tube cannot reach the stellar photosphere, but will be 
decelerated in situ, producing non-thermal (hard) X-ray emis- 
sion and collisional ionisation of elements such as Fe, which 
will be largely neutral at a typical temperature of ~ 7000 K 
dMuzerolle et all feOOll) . A sketch of the Elias 29 star-disk sys- 
tem with the hypothesised accretion tube where the electron ac- 
celeration is taking place is shown in Fig. [8] 

To explain the observed long-lasting fluorescent emission, 
the electron acceleration mechanism should be sustained over 
several days. While we have no detailed physical mechanism to 
propose for this, w e note that many of th e long-lasting flares ob- 
served in YSOs by Fav ata etaT] d2005al) require sustained heat- 
ing of the flaring material to explain the observed duration of the 
events, which in some cases lasted a few days. The thermody- 
namic cooling time of the flaring plasma would be much shorter; 
therefore, to explain the observed long flaring events in YSOs, 
magnetic reconnection must be ongoing, on time scales of days, 
similar to what is observed for the Fe 6.4 keV fluorescent emis- 
sion from Elias 29. The same electron acceleration mechanism 
could therefore result in long-lasting flares, if it takes place in 
systems with little ongoing accretion, or in strong Fe 6.4 keV 
emission, if it takes place in systems with strong ongoing accre- 
tion. 

In their study, Emslie et al. ( 1986) provide a relation (Fig. 1 
in their paper) between the flux emitted in the 6.4 keV line in the 
case of electron ionisation, as a function of the column density of 
material between the electrons' injection point and the "cold" Fe 
(where Fe 6.4 keV emission can take place). This column density 
is zero in our case, since, as mentioned above, the material in the 
accretion tubes is mostly neutral at temperatures of ~ 7000 K. In 
the diagram, different curves are given for different values of the 
spectral index, F, of the non-thermal X-ray emission, and the line 
flux is normalised to an amplitude of non-thermal (power-law) 
emission of a = 1 photons cm" 2 s" 1 keV" 1 (at 1 keV). 

The flux observed for Elias 29 in the Fe 6.4 keV emission 

line in the integrate d spectrum is ~ 10~ 6 photons cm" 2 s" 1 , 

and from the plot by lEmslie et all 1 1 9861 one can see that such 
a flux can be obtained for a very wide choice of parameter com- 
binations, with T ranging from 2.1 to 5.0, for a = 10" 4 - 1 
photons cm" 2 s" 1 keV" 1 . The parameters T and a refer to the 
hard X-ray em ission, typically at E <: 20 keV. As discussed by 
iHolmanl (12003). however, the X-ray spectrum of the non-thermal 
X-ray emission flattens below the low-energy cutoff (E c ) of the 
suprathermal electron responsible for the emission. For exam- 
ple, the bremsstrahlung emission from a beam of electrons, with 
a power-law spectrum with spectral index 5 = 3 and low-energy 
cutoff E c = 50 keV, has a spectral index of T ~ 2 at E > E c , but 
flattens to F ~ 1.3 between 1 and 10 keV. 

For the electron beam model in solar flares, electrons are 
assumed to have energy higher than 20 keV (e.g. iBrown et all 
1990) and recent RHESSI data indicates a low-energy cutoff 
that is typically E c = 20 - 40 keV and can be as high 70 keV 
dHolman et all 120031) . Using the BREMTHICK0 code by G. D. 
Holman, we verified that a non thermal hard X-ray emission with 
a = 10" 4 - 1 photons cm" 2 s" 1 keV" 1 and F = 2 - 5, from a pop- 
ulation of electrons with E c ^ 40 keV, between 1 and 10 keV, 
would be buried in the star's thermal emission and thus cannot 
be constrained without hard-X ray observations. Note, also, that 



3 available at http://hesperia.gsfc.nasa.gov/hessi/modelware.htm 
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Fig. 8. A cartoon of Elias 29 star-disk system with an accretion 
tube (hypothetical). A magnetic reconnection event within the 
accretion tube would generate accelerated electrons that could 
ionise the material and induce Fe Ka emission from the accret- 
ing material within the tube itself. The system is assumed at a 
viewing inclination angle i = 30° (i = 0° being face-on). 



a value of a — 10~ 4 - 1 photons ctrT 2 s _1 keV" 1 would not re- 
quire an unusually high power input in accelerated electrons. 

If our proposed mechanism of in situ electron bombardment 
of neutral or nearly-neutral material is correct, our observations 
provide additional evidence of strong non-thermal phenomena in 
young stars and, in particular, of particle bombardment of neu- 
tral material in the protoplanetary disk of the young star. In the 
solar system, such evidence is, for example, provided by the me- 
teorites known as chondrites. Chondrites have a number of pecu- 
liar features, as they contain granules (the chondrules) that have 
been flash heated (on time scales of less than one hour) to some 
2000 K, embedded in a rocky matrix with no evidence of heat- 
ing. Also, the chondrules show evidence of having cooled in a 
relatively intense magnetic field (of the order of 10 G). Finally, 
they contain material with peculiar isotopic ratios, with signif- 
icant amounts of short-lived radio nuclides (e.g. 41 Ca, with a 
half life of only 0.1 Myr, or 10 Be, with a half life of 1.5 Myr), 
which therefore cannot be originating in SN nucleosynthesis but 
must instead be originat ing in situ (as discussed in detail by 
e.g. lGounelle et al., 2006). While the observed isotopic anoma- 
lies likely require proton (rather than electron) bombardment, 
our observations provide evidence of in situ acceleration mecha- 
nisms closely associated with neutral material from the star's cir- 
cumstellar disk, that is, an environment where accelerated non- 
thermal particles interact with the protoplanetary material from 
which the meteorites will later grow. 

6. Conclusions 

So far, the Fe 6.4 keV emission in YSOs has been explained 
in terms of fluorescent emission from the photoionised (colder) 
material in the circumstellar disk. This scenario, however, cannot 
easily explain the observed variability of the Fe 6.4 keV emis- 
sion in Elias 29, which occurs in the absence of signficant varia- 
tions of the observed X-ray continuum. The equivalent width of 
the line at its maximum strength of ~ 250 eV is also not easily 
reconciled with a photoionisation scenario. An alternative line- 
formation mechanism is collisional excitation by a population 
of non-thermal electrons. We suggest that these electrons could 



be accelerated by magnetic reconnection events in the accretion 
tubes that connect the star to its circumstellar disk. The electrons 
are decelerated in situ by the accreting material, ionising it and 
causing the observed Fe 6.4 keV emission. 
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Table 1. Best-fit parameters from the fits to the spectra of Elias 29, during 6 different time intervals (PN data only), using an 
absorbed IT thermal model. Reduced^- 2 and null-hypothesis probability (P) for the fits are also given. 



Time interval 


AT(H) 


kT 


Z 


EM 


X 2 /d.o.f. 


P 


F x 




N 22 


keV 


Zo 


EM S3 






F-,3 


segl 


5.7 ± 0.7 


4.6 ± 0.9 


1.0 ± 0.3 


1.4 ±0.5 


1.0 


0.6 


1.8 


flare 


9.5 ± 2.0 


4.1 ± 1.5 


0.9 ± 0.3 


8.3 ± 6.3 


0.9 


0.5 


7.5 


seg2 


8.1 ± 1.3 


3.1 ±0.8 


1.0 ± 0.3 


1.4 ± 0.9 


1.5 


0.0 


1.1 


seg3 


6.6 ± 0.8 


3.4 ± 0.7 


0.7 ± 0.2 


1.5 ±0.8 


1.0 


0.5 


1.3 


seg4 


7.7 ± 0.6 


3.7 ± 0.5 


0.7 ±0.1 


3.1 ±0.9 


1.4 


0.0 


2.7 


seg5 


4.9 ± 0.4 


4.7 ± 0.8 


0.6 ±0.1 


2.0 ± 0.6 


1.3 


0.1 


2.5 



Units are N22 = 10 22 cm 2 , EM 53 = 10 53 cm 3 , and F_i 3 = 10 13 erg cm 2 s '. The spectral fits were carried out in the energy range 1.0-9.0 keV. 
Flux values refer to energy range 1.0-7.5 keV. The error bars are at lcr. 

Table 2. Best-fit parameters from the joint fits to PN, MOS1 and, MOS2 spectra of Elias 29, during 5 different time intervals, using 
an absorbed IT thermal model. 



Time interval 


2V(H) 


kT 


Z 


EM 


* 2 /d.o.f. 


P 


Fx 




N22 


keV 


z 


EM S3 






F-u 


seglf 


5.5 ±0.5 


4.4 ± 0.7 


1.0 ±0.2 


1.7 ±0.5 


1.1 


0.1 


2.1 


seg2t 


7.4 ± 0.8 


3.8 ± 0.7 


0.8 ± 0.2 


1.5 ±0.6 


1.6 


0.0 


1.3 


seg3 


6.6 ± 0.5 


3.5 ± 0.5 


0.7 ±0.1 


2.0 ± 0.5 


1.1 


0.2 


1.7 


seg4 


7.5 ± 0.4 


3.7 ± 0.4 


0.7 ±0.1 


3.5 ±0.8 


1.5 


0.0 


3.0 


seg5 


5.7 ± 0.4 


4.0 ± 0.5 


0.5 ±0.1 


2.8 ± 0.7 


1.6 


0.0 


2.3 



tMOSl data (mostly) unavailable for this time interval so only PN and MOS2 data were used. 

Table 3. Best-fit parameters from the fits to the spectra of Elias 29, during 6 different time intervals (PN data only), using an 
absorbed IT thermal model with an additional line at 6.4 keV (Fe fluorescent line). For the line, the total flux in the line (/6.4keV in 
units of /_6 = 10~ 6 photons ctrT 2 s _1 ) and the equivalent width (W(,a kev) are given. 



Time interval 


iVYH) 


kT 


Z 


EM 


/6.4keV 


^6.4 keV 


^ 2 /d.o.f. 


P 


Fx 


^6.4kevt 




N 22 


keV 


z Q 


EM 53 


f-6 


eV 






F-13 


% 


segl 


5.5 ± 0.7 


4.8 ± 1.0 


1.1 ±0.3 


1.3 ±0.5 


0.1 ±0.7 


13.0 


1.0 


0.5 


1.8 


45 


flare 


9.4 ± 2.0 


3.9 ± 1.4 


0.9 ± 0.3 


8.3 ± 6.2 


2.0 ± 5.9 


60.2 


0.9 


0.6 


7.5 


23 


seg2 


8.2 ± 1.3 


2.8 ± 0.6 


1.1 ±0.3 


1.4 ± 0.8 


1.1 ±0.6 


249.0 


1.3 


0.1 


1.1 


0.1 


seg3 


6.6 ± 0.8 


3.3 ± 0.6 


0.7 ± 0.2 


1.6 ± 0.8 


0.5 ± 0.5 


108.0 


0.9 


0.6 


1.3 


4.2 


seg4 


7.6 ± 0.6 


3.6 ± 0.4 


0.7 ±0.1 


3.2 ± 0.9 


0.7 ± 0.7 


67.8 


1.3 


0.1 


2.7 


3.5 


seg5 


4.9 ± 0.4 


4.5 ± 0.8 


0.6 ±0.1 


2.1 ± 0.6 


1.1 ±0.8 


162.0 


1.2 


0.2 


2.5 


0.2 


seg3+seg4+seg5« 


6.4 ± 0.3 


3.9 ± 0.3 


0.6 ±0.1 


2.2 ± 0.4 


1.0 ±0.3 


147.0 


1.3 


0.0 


2.1 


0.0 



^Probability that a line with a best-fit value equal or greater than We a kev could be due to a random fluctuation in the data, as estimated via Monte 
Carlo simulations (see Sect. 3 and 4). 

Results for the fit to the source's spectrum integrated over time intervals "seg3", "seg4", and "seg5" 

Table 4. Best-fit parameters from joint fits to PN, MOS 1 and MOS2 spectra of Elias 29, during 5 different time intervals, using an 
absorbed IT thermal model with an additional line at 6.4 keV (Fe fluorescent line). 



Time interval 


Af(H) 


kT 


Z 


EM 


/6.4keV 


Ws.4 keV 


* 2 /d.o.f. 


P 


Fx 




N22 


keV 


z 


EM 53 


f-6 


eV 






F_13 


seglf 


5.5 ±0.5 


4.3 ± 0.8 


1.0 ± 0.2 


1.7 ±0.5 


0.1 ± 0.7 


9.7 


1.2 


0.1 


2.1 


seg2t 


7.5 ±0.8 


3.5 ± 0.7 


0.9 ± 0.2 


1.6 ±0.6 


1.1 ± 0.7 


194.0 


1.4 


0.0 


1.3 


seg3 


6.7 ± 0.6 


3.3 ± 0.5 


0.7 ± 0.2 


2.1 ±0.7 


0.4 ± 0.6 


68.7 


1.1 


0.3 


1.7 


seg4 


7.5 ± 0.4 


3.5 ± 0.3 


0.7 ±0.1 


3.6 ±0.8 


1.2 ±0.7 


118.0 


1.4 


0.0 


3.0 


seg5 


5.8 ± 0.4 


3.7 ± 0.5 


0.5 ±0.1 


2.9 ± 0.8 


1.2 ±0.7 


152.0 


1.5 


0.0 


2.3 



tMOSl data (mostly) unavailable for this time interval so only PN and MOS2 data were used. 



